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■ >"H , 1 INTRODUCTION 



The continuum emission of stellar flares in UV and visible bands can be enhanced 
by two or even three orders of magnitude relative to the quiescent level and is usually 
characterised by a blue colour. Thermal atmospheric models are difficult to reproduce 
all these spectral features. If the flaring process involves the acceleration of energetic 
electrons which then precipitate downwards to heat the lower atmosphere, collisional 
excitation and ionisation of ambient hydrogen atoms by these nonthermal electrons 
could be important in powering the continuum emission. To explore such a possibility, 
we compute the continuum spectra from an atmospheric model for a dMe star, AD Leo, 
at its quiescent state, when considering the nonthermal effects by precipitating electron 
beams. The results show that if the electron beam has an energy flux large enough (for 
example, T\ ~ 10 12 ergcm~ 2 s _1 ), the U band brightening and, in particular, the U— B 
colour are roughly comparable with observed values for a typical large flare. Moreover, 
for electron beams with a moderate energy flux T\ ^ 10 11 ergcm~ 2 s _1 , a decrease of 
the emission at the Paschen continuum appears. This can explain at least partly the 
continuum dimming observed in some stellar flares. Adopting an atmospheric model 
for the flaring state can further raise the continuum flux but it yields a spectral colour 
incomparable with observations. This implies that the nonthermal effects may play 
the chief role in powering the continuum emission in some stellar flares. 
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Flares have been detected in stars of various types. In par- 
ticular, very energetic flares are likely to occur in dKe/dMe 
stars. The most energetic stellar flare can release an energy 
as large as ~ 10 35 erg in total, three orders of magnitude 
greater than that in the biggest flare on the Sun (see e.g. 
Pagano et al. 1997). Multiwavelength observing campaigns 
have been made to detect flare emissions in different wave- 
length regions like the soft X-ray, EUV, optical, and radio 
and then to quantify their relationships (e.g. Hawley & Pet- 
tersen 1991; Hawley et al. 1995; van den Oord et al. 1996). 
Though, there are still few hard X-ray observations of stellar 
flares. 

In spite of many different aspects between stellar and 
solar flares (Haisch, Strong & Rodono 1991), they may share 
similarities in some basic physical processes, such as the en- 
ergy release process, heating of the atmosphere, and origin 
of the continuum and line emissions. The well-known solar 
flare model assumes that a flare is produced through the re- 
connection of magnetic field lines, leading to the formation 
of a current sheet where tearing mode instability occurs. 



This process releases a large amount of thermal energy and 
results in an acceleration of electrons and/or protons. Accel- 
erated electrons then impact the lower atmosphere, produc- 
ing hard X-ray emission through bremsstrahlung, and caus- 
ing an enhancement of optical continuum and line emissions 
(e.g. Neidig et al. 1993) . Meanwhile, the corona accumulates 
material evaporated from the chromosphere when heated by 
nonthermal electrons and the soft X-ray flux gradually in- 
creases. The 'Neupert effect', assuming that the soft X-ray 
flux is proportional to the time integral of hard X-ray flux, is 
found valid not only for solar flares (Dennis & Zarro 1993), 
but also for stellar flares if using optical data as a proxy for 
hard X-ray emission (Hawley et al. 1995). 

The UV and optical continuum in stellar flares is found 
to exhibit several emission features. First, the continuum is 
usually characterised by a blue colour. Pagano et al. (1997) 
found in an intense flare a very blue emission (U — B w 
—2.5), implying a significant enhancement of the Balmer 
continuum and a flare site close to the stellar limb. A blue 
continuum was also observed by Hawley & Fisher (1992), 
which was reproduced from atmospheric models proposed 
by Mauas & Falchi (1996), taking into account the merging 
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of higher Balmer lines. Second, some flares are preceded by 
a dimming in their continuum emission (e.g. Hawley et al. 
1995). Several models have been proposed to explain such 
negative flares. Gurzadyan (1980) proposed inverse Comp- 
ton up-conversion of 'red' photons to the U band by ex- 
tremely energetic electrons. However, this model demands 
a total energy (~ 10 40 erg) of energetic electrons which is 
unrealistically large (Henoux et al. 1990). Mullan (1975) sug- 
gested that the dips in the red light curve can be produced 
by Ha absorption. At present, the most favourable model 
is the H~ absorption model of Grinin (1983), which sug- 
gests that an energy input to the lower atmosphere leads to 
an enhanced H~ opacity and then reduces the photospheric 
continuum radiation. 

The situation is more or less similar in solar white- 
light flares (WLFs), those flares emitting an enhanced opti- 
cal continuum. The majority of solar WLFs exhibit a pro- 
nounced Balmer jump while a small fraction shows no jump, 
which are classified as type I and type II WLFs, respectively 
(Machado et al. 1986; Fang & Ding 1995). Recombination 
of hydrogen atoms in the chromosphere is responsible for 
the continuum emission in the first case while the H~ emis- 
sion plays the main role in the second case. Besides, a neg- 
ative continuum contrast is also theoretically predicted at 
the very beginning of flare occurrence (Henoux et al. 1990; 
Ding & Fang 1996). Contrary to the stellar case, a negative 
solar WLF, the so-called black-light flare (BLF), is hard 
to observe due to its small amplitude of dimming and its 
short duration. In fact, no direct observational evidence has 
proved its existence (van Driel-Gesztelyi et al. 1994). 

Since the lower atmosphere of flares is very likely to 
be heated through bombardment of energetic electrons, the 
nonthermal collisional excitation and ionisation of hydrogen 
atoms can be important in increasing the continuum opac- 
ity and source function (Henoux, Fang & Gan 1993). The 
emergent continuum emission is then greatly enhanced and 
characterised by a Balmer jump. The continuum contrast 
depends on the flux and spectrum of the electron beam and 
also on the flare site (Ding & Fang 1996). 

The atmospheric structure and, in particular, the pa- 
rameters of the electron beam in stellar flares can be quite 
different from that in a solar flare. Thus, it is still unclear 
and needs to be quantified to what extent the nonthermal 
effects can alter the continuum emission in stellar flares. The 
purpose of this paper is to explore this question. To do so, 
we employ atmospheric models typical for dMe stars and do 
non-LTE calculations in the presence of an electron beam. 
The results show the role of nonthermal effects on the opti- 
cal continuum and are useful in the spectroscopy of stellar 
flares. 



2 COMPUTATIONAL METHOD 
2.1 Atmospheric models 

In last decades, attempts have been made to model the at- 
mosphere of stars of different spectral types. Here, we re- 
strict our attention to the work on late-type stars, in particu- 
lar, M dwarfs, which have a lower effective temperature. The 
reason is that in such stars, the afore-mentioned nonthermal 
effects during flares may be relatively more significant. 



Some early atmospheric models of dM or dMe stars were 
built on limited spectral data. Thus, their chromospheric 
structures are usually very schematic (e.g. Cram & Mullan 
1979). Multi-line spectral data are needed for a better mod- 
elling of the stellar atmosphere. Recently, after constructing 
grids of atmospheric models and examining their different 
spectral features, methods for probing the atmosphere are 
being developed which use spectral lines including H I Lya, 
Ha (H/3), Can H (K), Mgn h (k), and Nal D lines etc. as 
well as continua (e.g. Doyle et al. 1994; Houdebine & Doyle 
1994; Houdebine, Doyle & Koscielecki 1995; Houdebine et 
al. 1996; Houdebine & Stempels 1997; Short & Doyle 1998). 

Mauas & Falchi (1994) presented a semi-empirical 
model for the dMe star AD Leo, which can reproduce both 
the continuum and most spectral lines from observations. 
Their model is elaborate especially in the chromospheric 
structure, though, there may be some uncertainties in the 
temperature minimum region, as they stated. Mauas & 
Falchi (1996) further constructed semi-empirical models for 
a large flare on AD Leo observed by Hawley & Pettersen 
(1991). They computed three models corresponding to dif- 
ferent filling factors. Compared to the quiescent atmosphere, 
these models indicate a strong heating in both the chromo- 
sphere and the photosphere, which is quite different from the 
case of a solar flare. In particular, Mauas & Falchi (1996) 
discussed the possible role of high energy electrons on heat- 
ing the deeper atmosphere. 



2.2 Nonthermal excitation and ionisation rates of 
hydrogen by electron beams 

Here we assume that the basic scenario for solar flares is 
still ture for stellar flares, that is, flares originate in stellar 
coronae and involve the production of a beam of high en- 
ergy electrons, streaming downward to heat the lower atmo- 
sphere. In the cold target approximation, the rate of energy 
loss for beam electrons with column density N has the form 
(Emslie 1978; Henoux et al. 1993): 
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where K = 27re 4 . In equation (|l]), is the cosine of the 
pitch angle, £ is the ionisation degree of the target, A is the 
Coulomb logarithm due to collisions with target electrons 
and protons, while A' the one for inelastic collisions with 
target atoms. The latter corresponds to the energy deposit 
contributed to the nonthermal excitation and ionisation of 
target atoms. The rate of energy deposit at column density 
./V through inelastic collisions is written as (Emslie 1978; 
Chambe & Henoux 1979): 
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is the minimum energy of electrons that can penetrate to N. 
The cosine of the initial pitch angle, /io, is taken to be unity 
throughout the paper. The expression for /3 can be found in 
Emslie (1978). Knowing little about the energy distribution 
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of the beam electrons, we assume a power law for the initial 
flux spectrum with a low-energy cut-off: 



case A 



F (E) = 



AE~ 



E > Ei, 
E < Ei, 



(4) 



where A — (5 — 2)T\E 1 . T\ is the energy flux of electrons. 
Equation (@) then reduces to 
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is the column density penetrated by particles with an initial 
energy E\, and 



Mi = 



1, 

N/Nx, 



N > N u 
N < Nl 



(7) 



In the present calculations, we adopt an atomic model 
for hydrogen with four bound levels plus continuum. Ac- 
cording to Fang, Henoux & Gan (1993), the nonthermal col- 
lisional excitation rates due to bombarding electrons from 
the ground state are derived as 

C? 2 ~ 2.94 x 10 10 $/m, (8) 
d B 3 ~ 5.35 x 10 9 3>/ni, (9) 



and 

C? 4 ~ 1.91 x 10 9 <l>/ni, 

and the nonthermal ionisation rate is 

Cf c ~ 1.73 x 10 10 $/m. 



(10) 



(11) 



These rates are then incorporated into the statistical equi- 
librium equations, which are solved iteratively together with 
the radiative transfer equation for a specific atmospheric 
model. Finally the continuum flux at the stellar surface is 
computed as 

l 1 1 p 1 poo 

l'\ = 2- / I\(fi)fj,d^i — 2ty / &n 5 A e-^ /fl dr A . (12) 



3 RESULTS AND DISCUSSIONS 

Using the above method, we have studied the role of non- 
thermal electrons on the continuum emission in stellar flares. 
Since the flaring atmosphere undergoes a rapid evolution 
during the heating process, we consider two quasi-static 
cases (marked as A and B), an initial cool atmosphere and a 
heated atmosphere, which are assumed to be represented re- 
spectively by the atmospheric model for the quiescent state 
(Mauas & Falchi 1994) and that for a flaring state (Mauas & 
Falchi 1996, their model B) of the M dwarf AD Leo. This is 
reasonable as the radiative time scale is much shorter than 
the hydrodynamic time scale and our main interest is the 
radiative output. These two cases can be regarded as the 
very early and the maximum phases of flares respectively. 
The continuum emission is then computed for the two mod- 
els when considering the nonthermal collisional effects by 
electron beams precipitating from the corona. 
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Figure 1. Continuum spectra computed for a quiescent model 
(case A, left panel) and a flare model (case B, right panel) when 
considering the nonthermal collisional effects by precipitating 
electron beams. The energy fluxes, T\, of the beams are (solid 
line), 10 10 (dashed line), 10 11 (dotted line), and 10 12 (dash-dotted 



line) erg cm 



For all cases, the low-energy cut-off, E\ , is 20 



keV, and the power index, S, is 3. 

As an example, Fig. |l| plots the continuum spectra com- 
puted for both models in the presence of electron beams 
with a low-energy cut-off E\ — 20 keV and a power index 

5 = 3 but different energy fluxes T\ = 0, 10 10 , 10 11 , and 
10 12 ergcm _2 s _1 . It clearly shows that the bombardment 
of an electron beam leads to a dramatic increase of the con- 
tinuum emission in the Balmer continuum region (A < 3646 
A). The flux at the Paschen continuum region (A < 8204 A) 
changes to a relatively smaller extent. Therefore, a strong 
Balmer jump appears in the computed spectra which has 
not been observed (see also, Hawley & Fisher 1992). Mauas 

6 Falchi (1996) proposed that the merging of higher Balmer 
lines results in a pseudo-continuum which can smear out the 
Balmer jump. In the present case, the nonthermal electrons 
are expected to strengthen the Balmer lines as well (Fang et 
al. 1993). Thus, the merging effect could be even more im- 
portant here. In addition, the Doppler effect due to velocity 
fields of several hundreds of km s _1 possibly existing in the 
flaring atmosphere (e.g. Houdebine et al. 1993), and the line 
blanketing effect, which is not considered here, can also help 
to lessen the amplitude of the Balmer jump. 

3.1 Enhanced continuum emission and the 
spectral colour 

To make a comparison with observations, we extract the 
fluxes at A = 3600 and 4400 A, which lie in the centre of the 
Johnson U and B bands, respectively, and then compute the 
following parameters, 

AU' = -2.51og[i^3600,f/.F3600 J q], (13) 

AB' = -2.51og[F 44 oo,f/i 7 44oo, q ], (14) 
and 



A(U' - B') = AU' - AB'. 



(15) 



In the above equations, a subscript 'f refers to the flaring 
status (a quiescent/flare atmospheric model with a precip- 
itating electron beam), while 'q' to the pre-flare status (a 
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Figure 2. Computed parameters of At/', AB' , and A(U' — B') in 
dependence of the energy flux of the precipitating electron beam. 
For all cases, the low-energy cut-off, E\ , is 20 keV, and the power 
index, S, is 3. 

quiescent atmospheric model only). A prime means that the 
parameters are denned at fixed wavelength points instead of 
broad bands usually used in observations. Thus, the results 
can only be considered as exploratory. 

Fig. ^ displays the computed parameters versus the elec- 
tron beam fluxes. One finds that in the case of a large beam 
flux, pure nonthermal effects can help to raise the emission 
at A = 3600 A by two orders of magnitude. For example, if 
T\ reaches ~ 10 12 ergcm _2 s _1 , one obtains an enhance- 
ment corresponding to AU' ~ —5.3 (case A). Adopting 
a heated atmospheric model with the same electron beam 
yields further AU' ~ —6.9 (case B). This means that in the 
U band, the nonthermal effects could dominate over the ef- 
fect of temperature rise (the thermal effect) in producing 
the continuum emission. Pagano et al. (1997) reported on 
a large flare occurring on a dMe star G 102-21, which had 
a U band enhancement of « 7.3 mag. Taking into account 
the different pre-flare atmospheric conditions, we conclude 
that the U band emission in this flare can be qualitatively 
explained in terms of nonthermal effects. 

Under the current framework that the electron beam 
originates in the corona, it is conceivable that the Paschen 
continuum is less affected by the electron beam since it is 
formed in deeper layers than the Balmer continuum. Judging 
from Fig. the thermal contribution to the emission in the 
B band may remain to prevail over that due to nonthermal 
effects. For an electron beam with T\ ~ 10 12 erg cm -2 s , 
AB' is computed to be ~ —0.6 and —4.7 in cases A and B, 
respectively. The latter value is closer to the observations by 



Pagano et al. (1997), who obtained a B band enhancement 
of « 3.9 mag. 

Due to the greatly enhanced emission in the U band, 
the computed spectral colour appears very blue. The colour 
change, A (17' — B'), is computed to show the difference of 
relative enhancements at U' and B' wavelengths. In case A, 
this parameter varies in a broad range. The bluest colour, 

A(f7' - B') ' 4.8, appears when T\ ~ 10 11 erg cm" 2 s" 1 . 

The colour turns less blue if further increasing the values 
of T\. Thus, the observed colour change, A(J7 — B) « —3.4 
(Pagano et al. 1997), is easy to be reproduced. In case B, 
however, the spectra are contaminated by the thermal effect, 
the bluest colour of which is shown to be A(U' — B') ~ —2.3 
for electron beams considered here. The different results in 
these two cases imply that the observed blue colour is very 
likely of nonthermal origin. A pure thermal model is practi- 
cally hard to reproduce such a blue colour. 

Finally, it is worth noting that an electron beam can 
produce a return current instability when its flux exceeds 
some critical value. According to the formula derived by 
Aboudarham & Henoux (1986), we have found that electron 
beams with T\ ~ 10 12 erg cm -2 s _1 are stable in the corona 
of the flare model while they are only marginally stable in 
the corona of the quiescent model. Hence, a favourable sit- 
uation for the existence of such a large energy flux includes 
a relatively higher preflare coronal pressure and a harder 
spectrum of the electron beam (a smaller S and/or a larger 



3.2 The cause of a continuum dimming 

A continuum dimming is produced at the Paschen contin- 
uum (B' and V' wavelengths) in the case of a cool atmo- 
sphere (case A) bombarded by an electron beam with a mod- 
erate energy flux (T\ 10 11 erg cm -2 s" 1 ). However, such 
a darkening switchs quickly to a brightening with increasing 
values of T\. In the present models, we cannot produce a 
dimming at the Balmer continuum. 

The continuum dimming results from attenuation of the 
photospheric radiation by an enhanced opacity in the lower 
chromosphere and the temperature minimum region. Non- 
thermal ionisation of hydrogen atoms overpopulates the am- 
bient electrons, leading to an enhanced H~ opacity. There- 
fore, the magnitude of the depression depends sensitively 
on the electron energy, which determines the electron pen- 
etration depth. Fig. fa illustrates the effect on the contin- 
uum dimming at A = 4400 A of varying the low-energy 
cut-off, while keeping the energy flux of the electron beams 
unchanged. It clearly shows a positive dependence of the 
magnitude of the dimming on the electron energy. 

In observations, the U band comprises a part of the 
Paschen continuum. Thus, we think that the observed con- 
tinuum dimming could be related to an electron beam bom- 
barding a relatively cooler atmosphere. This is likely to occur 
in the early phase of flares. In conclusion, the nonthermal 
effects could at least partly, if not fully, account for the con- 
tinuum dimming observed in some stellar flares (e.g. Hawley 
et al. 1995). 
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Figure 3. Computed continuum dimming, -F4400 f/^4400 qj m 
dependence of the low-energy cut-off of the precipitating electron 
beam. For all cases, the energy flux, Ti, is 10 11 ergcm _2 s -1 , 
and the power index, 8, is 3. 



3.3 Interpretations of the results 

To understand more clearly the above results, we give in 
Fig. ^the height distributions of the source function, optical 
depth and the contribution function per unit geometrical 
depth, 



dJ\/dz = 2tt 



\./a» 



(16) 



at A = 3600, and 4400 A for case A. It shows that the 
nonthermal effects can raise both the source function and 
the opacity to a large extent in the chromosphere. Thus, 
the chromospheric contribution to the emergent flux gradu- 
ally increases while the photospheric contribution decreases 
with increasing electron beam flux. Whether the continuum 
brightens or darkens is then determined by the competition 
of these two factors. 



4 CONCLUSIONS 

If stellar flares result from the reconnection of magnetic 
fields, analogously to solar flares, electrons and/or protons 
could be accelerated at the reconnection site and then pre- 
cipitate into the lower atmosphere. The nonthermal colli- 
sional excitation and ionisation of ambient hydrogen atoms 
by these electrons lead to an enhancement of the source 
function and opacity in the atmosphere. The chromospheric 
contribution to the continuum emission rises rapidly, pro- 
ducing a spectrum with a very blue colour. By employing 
an atmospheric model for an M dwarf star, AD Leo, at 
quiescent state (Mauas & Falchi 1994), we have computed 
the continuum spectra for various electron beams bombard- 
ing the atmosphere from the corona. The results show that 
an intense electron beam with an energy flux, for example, 
J-\ ~ 10 12 ergcm _2 s _1 can produce a continuum flux at 
A = 3600 A (Balmer continuum) two orders of magnitude 
greater than the quiescent value. The U band brightening 
and, in particular, the U — B colour are roughly compara- 
ble with observed values for a typical large flare. Adopting 
an atmospheric model for the flaring state can further raise 
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Figure 4. Height distributions of the source function, optical 
depth and the contribution function at A = 3600, and 4400 A for 
case A. The energy fluxes, T\, of the electron beams are (solid 
line), 10 10 (dashed line), fO 11 (dotted line), and fO 12 (dash-dotted 
line) erg cm -2 s — 1 . For all cases, the low-energy cut-off, E\, is 20 
keV, and the power index, S, is 3. 



the continuum flux and can better account for the B band 
enhancement. However, it yields a spectral colour incompa- 
rable with observations. This implies that the nonthermal 
effects may play the chief role in powering the continuum 
emission in some stellar flares. 

Computations also predict a continuum dimming at the 
Paschen continuum, due to attenuation of the photospheric 
radiation by an enhanced opacity in the lower chromosphere 
and the temperature minimum region. This happens in a 
cool atmosphere bombarded by an electron beam with a 
moderate energy flux {T\ 10 11 ergcm _2 s _1 ). This is 
likely to fit the conditions in the very early phase of flares. 
The magnitude of the dimming depends positively on the 
electron energy. We argue that the continuum dips found 
in observations (e.g. Hawley et al. 1995) could possibly be 
related to nonthermal electrons of higher energies. 
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